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Abstract 

Accurate knowledge of the single-molecule (self) translational dynamics of liquid para-H 2 is an 
essential requirement for the calculation of the neutron scattering properties of this important 
quantum liquid. We show that, by using Centroid Molecular Dynamics (CMD) quantum simula¬ 
tions of the velocity autocorrelation function, calculations of the total neutron cross section (TCS) 
remarkably agree with experimental data at the thermal and epithermal incident neutron energies 
where para-H 2 dynamics is actually dominated by the self contributions. This result shows that 
a proper account of the quantum nature of the fluid, as provided by CMD, is a necessary and 
very effective condition to obtain the correct absolute-scale cross section values directly from first- 
principle computations of the double differential cross section, and without the need of introducing 
any empirically adjusted quantity. At subthermal incident energies, appropriate modeling of the 
para-H 2 intermolecular (distinct) dynamics also becomes crucial, but quantum simulations are not 
yet able to cope with it. Existing simple models which account for the distinct part provide an 
appropriate correction of self-only calculations and bring the computed results in reasonable accord 
with TCS experimental data available until very recently. However, if just published cross section 
measurements in the cold range are considered, the agreement turns out to be by far superior and 
very satisfactory. The possible origin of slight residual differences will be commented and suggest 
further computational and experimental efforts. Nonetheless, the ability to reproduce the total 
cross section in the wide range between 1 and 900 meV represents an encouraging and important 
validation step of the CMD method and of the present simple algorithm. 

PACS numbers: 61.05.F, 83.10.Rs, 67.63.Cd 
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I. INTRODUCTION 


Hydrogen isotopes H2 and D2 in their liquid states have always attracted much inter¬ 
est for their importance as fundamental and simplest molecular systems. The “hydrogens” 
display a quantum behavior still challenging theoretical and simulation-based descriptions 
of their translational dynamics, as far as both single-molecule and collective properties are 
concerned. The study of these liquids is, therefore, evidently important per se; nevertheless, 
its relevance is even further enhanced by two facts. One is the direct access to their micro¬ 
scopic dynamics enabled by important and widespread spectroscopic techniques. Neutron 
scattering, in particular, is by far the most important means for studying the hydrogens and 
hydrogen-containing materials, and it is well known that the study of the cross section for 
scattering of neutrons from a sample of liquid hydrogen has been, by itself, a long-standing 
theme in the physics of quantum fluids (see e.g. [ 1 ] and Refs, therein). The second fact is 
that liquid hydrogens are among the most used cryogenic fluids and, in the specihc appli¬ 
cation to neutron techniques, are the most important low-temperature neutron moderators 
used to realize cold neutron sources. 

In this field of application, an improved description of the neutron double differential cross 
section (DDCS) and total cross section (TCS) of liquid H2 and D2, and of other moderating 
materials, has become an indispensable requirement for appropriate development and up¬ 
grade of neutron facilities. For instance, it has been recently showrPthat calculations of the 
safety rod insertion impact on the criticality of water-moderated reactors provide much bet¬ 
ter results if based on new accurate (experimental and/or simulated) DDCS determinations, 
rather than on existing cross section libraries and scattering kernels employed in the available 
nuclear data processing codes treating water and other neutron moderators.®^ Similarly, the 
capability of predicting, with a high reliability, the dynamic response to neutrons of the hy¬ 
drogen liquids is of crucial importance for the design of advanced and high brilliance cold 
neutron sources, as those aimed to exploit directional moderator geometries.^ 

Direct experimental investigations of the scattering law of these systems are of course 
necessary to validate any model calculation, but obviously cannot cope with the present need 
to know the DDCS in a large number of conditions of exchanged wave vector Q and energy 
E for any relevant incident energy Eq. Therefore the development of computable accurate 
models for the DDCS is a major objective in this held of research and application, especially 


3 


if focused on first-principle methods aimed to obtain direct agreement with experimental 
data without forcing it through the introduction of variable parameters or ad hoc values 
for some physical quantities. Unfortunately, available data, appropriately normalized to an 
absolute scale, mainly refer to integrated quantities like the TCS measurements of Refs. [ 8 ^ 
[To], while past DDCS determinations for hydrogen and deuterium are often given in arbitrary 
units [see e.g. ITT1TT3] and exclude the possibility of a “complete” (i.e. both in shape and 
absolute-scale intensity of the spectra) comparison between data and calculations. 

The construction of appropriate models is however made rather complex by the quantum 
nature of the hydrogen liquids, which is highlighted by the relatively large values of the de 
Broglie thermal wavelength A = h/\/ 2 'KMk-QT, where M is the molecular mass, T is the 
temperature, and k-Q and h are the Boltzmann and Planck constants, respectively. Indeed 
the low mass and low temperature of these liquids makes A comparable with the molecu¬ 
lar dimensions, so that distinguishability of the particles is retained (and use of quantum 
statistics is unnecessary) but delocalization effects show up in the scattering properties, as 
evidenced even for the heavier deuterium molecule (see e.g. dl)- 

Therefore classical analytical or simulation methods are not appropriate to predict the 
center-of-mass (CM) single-molecule {self) and total {self plus distinct) dynamic structure 
factors of liquid H2 and D2. It is important to stress that the hydrogens are “quantum” 
liquids in the above special meaning, which is weaker than that of the He case, but is anyway 
more profound than that brought about by the non-commutative properties of position and 
momentum operators, leading to detailed balance spectral asymmetry and non-zero hrst 
frequency moment even in a non-interacting monatomic system. Furthermore, quantization 
of the internal degrees of freedom (e.g. rotations and vibrations) of molecular fluids is once 
again something else, which has to be duly considered also for “classical” liquids such as 
water or methane, and not only in the case of the hydrogens. 

Much work was devoted in past years to the implementation of scattering kernels for 
the hydrogen liquidP^^^^ with different methods used to effectively evaluate both the self 
and distinct CM contributions to the DDCS. However, the above mentioned quantum effects 
were mostly neglected, or accounted for only approximately, or hnally included in an effective 
way by adjustments to pioneering experimental data.^nHl 3 | i 9 | 20 | 

In order to investigate their role in the neutron response of these liquids, we focus here 
on the most signihcant case of H2 at liquid temperatures (i.e., para-H2), which combines a 
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stronger quantum behavior with a full predominance (in most kinematic conditions) of the 
self component, the latter being quite reliably accessed by quantum simulation methods (see 
Sect. Ill) and sensitively probed by neutrons in a unique way. 


II. THE DDCS OF Ha 


By assuming free rotations and neglecting rotation-vibration coupling, the nuclear neu¬ 
tron scattering by a homonuclear diatomic molecule can be schematically described by 


withffl 


dVtdE 



Sn{Q,E) 


S^{Q,E) = u{Q)ScmMQ,E)+ 

+ ^ Fjgj-^y-^{Q)ScM,seifiQ^ E — Ejgj-^^ — Eo^-^) ^ ^ 

JqJivi 

where u{Q) is a Q-dependent function containing only the coherent cross sections of the 
nuclei in the molecule, and the function E takes different expressions according to the 
nuclear spin statistics and the ortho-para concentration, and contains both the coherent 
and incoherent nuclear cross sections (see H] ). In Eq. (|5 ^CM,dist(Q,^) and ^cM,seif(Q, ^) 
denote the distinct (j 7^ j') and self (j = j') components of the total dynamic structure 
factor per molecule Scm{Q, E): 
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( 2 ) 


ij'=i 


where N is the number of molecules and Rj and Rj/ denote the CM position operators of 
the j-th and /-th molecule at time 0 and t, respectively. 

The second term of Eq. ([^ represents the single-molecule dynamics S'cM,seif(Q, T^) con¬ 
voluted with the line structure of the internal molecular motions, and results in a sum of 
spectral lines centered at the energies of rotational (Jq —t Ji) and vibrational (uq —t Vi) 
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transitions, where the snbscripts 0 and 1 are nsed to label initial and hnal state, respec¬ 
tively. The gronnd vibrational state (wq = 0 ) is assnmed here as the only one signihcantly 
popnlated in hydrogen at liqnid temperatnres. 

As anticipated, the last term in Eq. Q is particularly important in the case of liq¬ 
uid H2 due to the huge incoherent-to-coherent cross section ratio of the hydrogen nucleus. 
Moreover, in the case of para-H2 the combination of initial-state probabilities and spin cor¬ 
relations makes the DDCS spectra dominated, in most kinematic conditions, by the intense 
rotational line with Jq = 0 and Ji = 1 , at ~ 14.7 meV.^^ Consequently, a good model¬ 
ing of S'cM,seif(<5, .E) alone is generally sufficient to account for the scattering of thermal 
and hot neutrons from this liquid. In this respect, the results on liquid para-H2 provide a 
fundamental test of the general quality of possible representations for the single-molecule 
dynamics. 

The simplest algorithms for the evaluation of the single-molecule part of the 
originate from direct or empirically modified use of, basically, two possible analytical models 
for S'cM,seif(Q, either the ideal gas (IG) law (see e.g. [ 1 ]), giving rise to the well-known 
Young and Koppel modeP^for the self DDCS; or the Egelstaff and Schofield (ES) model,!^ 
modified to comply with detailed-balance asymmetry,!^ and with the first frequency moment 
sum rule that ensures translational spectra with a (non-zero) first moment equal to the recoil 
energy / 2 M. By contrast, in Ref. [ 15 ] use is made of a semiempirical determi¬ 

nation of the spectrum of the CM velocity autocorrelation function (VACF) of a hydrogen 
molecule, which, in the Gaussian approximation (GA),I 2 ZE 3 jg related to the intermediate 
scattering function TcM,seif(QC) and, through a time Fourier transform, to S'cM,seif(Q,-E). 
However, in Ref. [TS| the GA was not applied in the rigorous quantum-mechanical form 
given, for example, in Ref. [ 2 B] and summarized in the next section. 

Although different from each other, all the above approaches suffer anyway from the 
limitation of neglecting or approximating the quantum behavior of the liquid, in the sense 
we specified previously. Moreover, none of them provides direct agreement with experiment 
unless by modifying some property entering the model (like e.g., in the ES model, the mass 
of the molecule or the self diffusion coefficient) or by resorting to adjustable parameters in 
the spectrum of the VACF, as in the case of Ref. na). In some cases, TCS measurements 
have been fairly reproduced by considering the effect of intermolecular vibrations on the 
VACF spectrum and by hypothesizing the presence of huge molecular clusters diffusing in 
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liquid hydrogen, leading to a ~ 35-40 times augmented molecular to be used 

in the diffusive part (ES) of the single-molecule dynamics. However, such a high mass value, 
usually explained in terms of translational hindering and cluster formation, is not easily 
justihed, and might reveal instead some inadequacy of the ES lines-hape. 

These observations induced us to explore whether model parametrization and fit-based 
adjustments to experimental data could be fully avoided by means of a quantum determi¬ 
nation of S'cM,seif(Q, E). Moreover, this could help solving known physical inconsistencies, 
like the violation of the quantum second moment sum ruld® by the ES and IG spectral 
line-shapes. In addition to the compliance with the quantum properties of the liquid, other 
advantages of such a parameter-free and first-principle method are evident: i) a straightfor¬ 
ward adaptation to different thermodynamic and kinematic conditions; ii) the possibility to 
avoid introducing not well-justified hypotheses and to use the true H 2 molecular properties; 
iii) good control of the physical consistency of the used line-shape for S'cM,seif(Q, E) and of 
its compliance with the basic sum rules for a quantum fluid. 

The next section illustrates the method and the successful achievement of the above 
expectations both for the spectral properties of 5'cM,seif(Q, E) and for the neutron TCS 
results at thermal and epithermal neutron energies. 


III. QUANTUM SIMULATION-BASED SELF DYNAMICS 

The rationale behind the attempt here described is given by the recent validation of 
quantum Centroid Molecular Dynamics (CMD) as an effective method to simulate the VACF 
of hydrogen and the detailed probing of the degree of accuracy (and range of applicability) 
of the Gaussian approximation in predicting the measured line-shape of liquid para-H 2 .^ 
Indeed, a non-Gaussian behavior was observed in an intermediate Q range. However, the 
differences between neutron data and calculations, combining CMD simulations with the 
GA (denoted in the following as CMD-I-GA), are small enough to suggest that they may be 
irrelevant for cross-section calculation purposes, especially at the level of doubly integrated 
quantities like the TCS. 

As described in Ref. EH, the Path Integral CMD method was applied to a system 
of 256 molecules interacting via the Silvera-Goldman potentiaP^l. The Trotter number, 
i.e. the number of beads on the classical ring polymers replacing the quantum mechanical 
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particles, was 64. In contrast to the usual implementation, the calculation of the quantum 
mechanical forces, which are required at each time step of the otherwise classical simulation, 
was performed by the path integral Monte Carlo method, rather than MD, thus avoiding 
sampling problems associated with the stiff “intramolecular modes of the polymers and 
allowing for a much larger time step. The simulation was extended up to 1 ns in the isokinetic 
ensemble, ensuring thermal stability and statistical reliability. The velocity correlation was 
calculated up to a maximum time lag of 1.5 ps. A shorter test run with 500 particles 
conhrmed that the shape of the VACF was not noticeably influenced by hnite-size effects. 

The dynamical information conveyed by the VACF is a keypoint in the development of 
models for the self part of the DDCS of viscous dense fluids. In particular, the output of a 
PI CMD simulation is the canonical (or Kubo-transformecP^ VACF: 


Uc{t) = k^T 


dX 




vcm(O) 




Vcm(^) 


( 3 ) 


Jo \ j 

where H is the Hamiltonian operator of the system. The self intermediate scattering function 
-^CM,seif(Q, in the Gaussian approximation can then be written a^^ 
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Figure[^reports our simulation results for the VACF of para-H 2 at 15.7 K: both the canonical 
VACF Uc(t) and the real part of {vcm(O) ■ vcM(t)) are displayed for comparison with the 
VACF values digitalized from Fig. 3 of Ref. |T5], also shown in the hgure. The differences 
between the present CMD results and the past semi-empirical determination of Ref. [15] at 
14.7 K seem too large to be ascribed to a mere temperature effect. 








FIG. 1. (Color online) Velocity autocorrelation function of hydrogen deduced in [15] (pink dots). 
The quite different results of the quantum CMD simulations^ are also reported. For comparison, 
both the canonical VACF Uc{t) (dash-dotted blue curve) and the real part of (vcm(O) • ’vcuit)) 
(red solid curve) are shown. 

The self dynamic structure factor, as a function of uj = E/h, is then obtained as the time 
Fourier transform of Eq. (|^ at each desired Q. By integration of the spectra over very 
wide energy ranges, it is possible to verify the consistency with the second moment sum rule 
which, for a quantum system, reads® 


M(2)(Q) = JdE 5cM,self(Q,E) = 


3M 


{Ek) + E^, 


( 5 ) 


where (Ek) is the mean kinetic energy of the particle, which in the present case differs 
signihcantly from the classical value (3/2) k^T. Experimental and Path Integral Monte 
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Carlo simulation values of (-Er) for para-hydrogen at various liquid temperatures have been 
provided by Celli et al® and Colognesi et ah.l^ Fig. shows the Q-dependence of the 
theoretical prescription of Eq. (|^ (lines) in comparison with the IG, ES and CMD-I-GA 
results for the second frequency spectral moment (symbols). In particular, the continuous 
red curve corresponds to the calculation of Eq. ([^ using the experimental value of the mean 
kinetic energy of para-H 2 at 15.7 K,^ while the dashed black curve is derived by assuming 
a classical mean kinetic energy of (3/2) k^T. As expected, the latter calculation agrees very 
well with the values obtained by appropriate energy integration of the IG spectra. Similarly, 
the GMD-I-GA values are in very good agreement with the quantum behavior. Gonversely, 
the ES model misses both the classical and quantum prescriptions of 


IV. COMPARISON WITH EXPERIMENTAL DATA IN THE THERMAL AND 
EPITHERMAL RANGE 

The previous results show that the GMD-I-GA S'cM,seif(Q, E) has the correct spectral 
properties and can sensibly be used to perform the self DDGS estimates based on the com¬ 
putation of the second term of Eq. ([^, which we implemented according to the details given 
in Ref. [1]. By numerical integration, over energy and solid angle, of the DDGS spectra we 
hnally derived the TGS values reported in Fig. which, as mentioned, regard an incident 
energy range where the distinct contributions to Eq. ([^ are fully negligible to a good approx¬ 
imation. Galculations were carried out also using the IG and ES line-shapes for comparison 
with the quantum simulation-based outputs. From the inset of Fig. the superiority of 
the GMD-I-GA results in reproducing the experimental TGS of para-H 2 at thermal incident 
energies is evident. At higher energies, as expected, the system tends to ideal gas behavior: 
all calculated curves thus become indistinguishable and in very good agreement with the 
experimental data. 

At the more detailed level of non-integrated quantities, the effectiveness of the GMD-I-GA 
method can be partially tested against one of the few inelastic scattering measurement^ 
providing data in absolute units. In this case, since the hydrogen sample had a 41% ortho 
concentration, the GMD-I-GA computations were performed in the same conditions. The 
GMD-I-GA calculated spectra satisfactorily describe the experimental data of Fig. This 
shows that use of a quantum representation of S'cM,seif(Q, E) is the only real requirement 
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FIG. 2. (Color online) Second frequency moment of para-H 2 at 15.7 K calculated from Eq. ([^ 
either using the experimental (Ek) estimate of Ref. |34j (red solid curve), or the classical mean 
kinetic energy (3/2) k-^T (black dashed curve). Symbols correspond to the values obtained instead 
by energy integration of the IG (pink diamonds), the ES (blue squares), and the CMD+GA (black 
dots) spectra. 

to obtain a direct and reasonable agreement with experiment, comparable with the one 
obtained with the calcnlations of Ref. [TS], and incomparably better than that of Schott 
himself (using the IG model) and of UtsurcP^ (both not shown in the hgure). 

The method of combining VACF quantum simulations with the GA represents therefore 
an extremely valid alternative to most experiments or to the use of “quantum-insensitive” 
analytical models for S'cM,seif(Q, .F) in the DDGS formula of H 2 . A great advantage of this 
method, recently applied with simpler classical simulations to the case of wateris that the 
VAGF only depends on the thermodynamic state, and only one simulation run is required 
to enable, in a given state, calculations of S'cM,seif(Q,-E) wherever wished in the kinematic 
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FIG. 3. (Color online) Total cross section per molecule of para-H 2 at 15.7 K and 1 atm, at the 
thermal energies of Seiffert measurement^ (black dots in the inset) and at the higher Eq values 
probed in the spallation-source experiment of Celli et al.® (black empty circles in the main frame). 
Experimental data are compared with the IG (pink dash-dotted line), ES (blue dashed line) and 
CMD-I-GA (red solid line) calculations. 


{Q,E) plane. This actually means that, at least for H 2 , we are presently able to calculate the 
DDCS in most kinematic conditions, with good accuracy, and high control on the physical 
consistency of the results. 


V. COMPARISON WITH EXPERIMENTAL DATA IN THE SUBTHERMAL 
RANGE 

Differently from the case of more energetic incident neutrons, at cold neutron energies 
(1 < < 10 meV) the interplay of initial-state probabilities and spin correlations in para- 

hydrogen actually kills out the incoherent signal, leaving the small coherent cross section at 
the rudder of both collective and single-molecule contributions. Indeed, if only quasielastic 
scattering is allowed (because cold neutrons are unable to induce transitions in para-H 2 ) 
and only the Jq = 0 level is thermally populated (because of the low temperature of this 
liquid) it can be shown that the second term of Eq. ([^ reduces to u{Q)ScM,seuiQ^ that 

S'n(Q, E) = m(Q)S'cm(Q, E). As a result, at neutron energies below the threshold of the hrst 
rotational transition for this system, the distinct dynamics warrants a primary role also in 
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DDCS [bam/sr/meV/molecule] 



Scattered neutron energy Ei [meV] 


FIG. 4. (Color online) Experimental DDCS of H 2 at 19.8 K with a 41% ortho concentration, (black 
squares)!^ and CMD+GA results (red dots with thin line). The blue empty circles are the DDCS 
values calculated by Morishima and Nishikawa.^^ Measurements refer to the thermal range {Eq = 
21.8 meV) and rather wide angles indicated in each frame. 


the response from a nominally ’’incoherent” liquid as para-H 2 . Actually, there are evidences 
that distinct intermolecular contributions influence the neutron signal from para-H 2 and 
that self calculations are insufficient in certain conditions.^ 

A good knowledge of the collective (i.e. coherent) dynamics of this liquid is therefore 
an essential ingredient to achieve the accuracy demanded nowadays on subthermal neutron 
cross sections. Unfortunately, present quantum simulations methods are not yet able to 
provide direct and reliable estimates of the self and total (self plus distinct) dynamic struc¬ 
ture factors.l^^® At the same time, no analytical model exists for the total ScuiQiE) of a 
quantum liquid, and more generally of any liquid, except that in hydrodynamic {Q —)■ 0) 
conditions. Therefore, we could only try to investigate, in an approximate way, the effect 
of adding a collective term to the single-molecule results. To do this, we adopted the Skold 
approximation,^ that models the total ScuiQ, E) through a modihcation of its self part, 
namely: 


^cm(Q,^)~^cm(Q)^cm ,self 


Q 

\/ Scm{Q) 



( 6 ) 


which was cleverly conceived to fulhll the second moment sum rule in classical systems, 
and gave satisfactory results in several cases, starting from that of liquid argon.®^ In Eq. 
(§ Scm{Q) is the CM static structure factor for which we took the experimental values 
obtained by neutron diffraction measurements.®^ 
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As evident in Fig. the CMD+GA calculation of the self part progressively departs 
from the TCS measurements presently availabl#*^, as soon as the incident neutron energy 
is decreased below 14 meV. The Skold schematization of the total dynamics, calculated by 
inserting in Eq. (§ the properly Q-scaled CMD+GA self dynamic structure factor, is seen 
instead to be rather effective, providing a suitable correction of the self-only results. In 
particular, despite its simplicity, the Skold model, if used in combination with a quantum 
*S'cM,seif(Q 5 -E), is able to quite satisfactorily describe the very recent neutron transmission 
data of Grammer and co-workers at 15.7 K, with only a small overestimate between 2 and 
8 meV. Our quantum-based calculations thus apparently conhrm the lower cross section 
values of liquid para-H 2 with respect to the earlier cold neutron determinations by Seiffert. A 
possible contamination, by a small fraction (0.5 %) of ortho molecules, of Seiffert sample has 
been plausibly hypothesized in Ref. [10]. This actually seems to be the case, as shown in Fig. 
1^ where we also report, for an example incident energy of 3 meV, the changes induced, in 
both the GMD+GA and Skold TGS, by the presence of 0.5 % ortho-H 2 . Indeed, temperature 
effects might also be thought of to hnd an explanation to the quite large discrepancies 
between the two TGS data-sets (15.7 K against 14 K). However, it is well known that 
structural changes are extremely limited with varying temperature in high density liquids,®' 
and could unlikely lead to such signihcant differences in the TGS. 

The origin of the residual differences between Grammer et al. data and our calculations 
in the range 2< Eq <8 meV, cannot be easily found out. They might either depend on 
inaccuracies of the data, or of the S{Q) determination at the Q values involved at certain 
incident energies, or hnally on both. It is also plausible that the featureless line-shapes 
produced by the Skold recipe do not provide an adequate description of the total dynamic 
structure factor which, besides the single-molecule dynamics, accounts for collective excita¬ 
tions as well, with well-known inelastic features. At low and intermediate wave vectors (e.g. 
up to that of the main peak of S{Q)), substantial spectral differences from the Skold model 
are therefore expected, and might reveal themselves in a constant-^ energy integration be¬ 
tween (ideally) —oo and Eg, like the one performed to compute the TGS. Finally, another 
possible explanation of the small deviations could be a slight inadequacy of the GMD+GA 
method, and in particular a signature of the mentioned non-Gaussian behavior of para-H 2 
in certain Q ranges. In this respect, some observations are worth. Although deviations from 
Gaussian behavior were found to be rather small at the level of frequency spectra,®! the 
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FIG. 5. (Color online) Total scattering cross section of para-H 2 at cold neutron energies. The 
large deviation of the CMD+GA self contribution (red solid curve) from the absorption-corrected 
experimental data of SeifferiP (black dots) and of Grammer et (black empty circles) is due to 
the missing negative distinct contribution. The pink dotted curve displays the results obtained by 
adding a distinct part in the Skold approximation. The dots a 3 meV show the effect, on the self 
(red dot) and on the self-l-Skold (pink square) TCS, of the presence of 0.5 % ortho-hydrogen. 


effect of even small discrepancies has a different impact on the TCS results depending on 
and is more important if the Q values where the GA fails belong, in majority, to the 
explored kinematic region. Conversely, at very low Eq (e.g. 1-2 meV) the probed Q region 
predominantly lies between values where the GA holds, explaining the renewed agreement, 
as energy is decreased below 2 meV. Finally, at higher Eq (e.g. between 8 and 14 meV) 
many, rather low and rather high, Q values (where the GA holds) contribute to the TCS, so 
the role of the “GA-deviating” Q values has a much lesser impact on the integrated result, 
and agreement with experiment is found again. In our opinion, this might be a sensible 
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explanation of the initially growing and then decreasing discrepancies with increasing Eq in 
the snbthermal range. Anyway, it is notefnl that the present CMD+GA approach allows to 
reach, nsing in the calcnlations the true molecular properties of para-H 2 , such a satisfactory 
description of the new measured data. 

Clearly, these are plausible speculations that can only be verihed by further research. 
Two different routes can be envisaged in order to understand, if possible, the role of the 
two main approximations here adopted (i.e., GA and Skold) to calculate the snbthermal 
TCS. One is the introduction in our DDCS algorithm of the hrst non-Gaussian correction 
to the self dynamics.®* Though feasible, such an attempt is certainly “expensive” from the 
computational point of view, but it is also extremely attractive in a scientihc sense, owing 
to the great opportunity to test possible (near to “macroscopic”) failures of the GA, in case 
they show up even at the level of neutron TCS calculations. In addition, the dependence 
on Eq of the supposed GA inadequacy could be duly checked. The other route passes 
through experiments aimed at determining the dynamic structure factor and comparing it 
with the Skold prediction. Indeed, accurate measurements of the DDCS of liquid para-H 2 , 
used in combination with the available simulations of the self part, are needed to be able to 
sensitively probe the different contributions to the total signal in the sub-thermal and cold 
incident energy range below ~ 14 meV, with a great effort in producing well-normalized and 
duly corrected scattering data. 


VI. CONCLUSIONS 

We showed that it is now possible to obtain accurate evaluations of the thermal neutron 
scattering law of cryogenic liquids as important as para-H 2 . We developed, verihed and im¬ 
plemented an efficient simulation-based method able to accurately account for the quantum 
behavior of the huid, and which has the doubtless merit of limiting considerably the need of 
consuming experiments, at least as regards the thermal and epithermal range. The whole 
kinematic plane can indeed be covered by this technique, avoiding demanding experiments 
and with full hexibility. First-principle methods, use of the known molecular properties of 
hydrogen, and neither search of forced agreement with experiment nor introduction of ef¬ 
fective mass values or other parametrized routes, have been successfully experimented. An 
improved evaluation of the non-Gaussian behavior of para-H 2 is suggested by the present 
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TCS results, which point also at the importance of accurate measurements of the DDCS at 
some subthermal incident neutron energies where the role of the distinct contribution is, as 
well, worth investigating. The present estimates represent, to our knowledge, the only ones 
providing a parameter-free absolute-scale agreement with total cross section measurements, 
ensuring, at the same time, compliance with the spectral quantum sum rules. In this sense, 
the results reported in this work also constitute a rather convincing validation, at present, 
of the CMD simulation technique for the prediction of the VACF of quantum liquids. This 
also highlights the importance of quantum simulation research, both for fundamental mo¬ 
tivations, and for its fruitable use in increasingly demanding helds of application, like the 
management and new-concept development of neutron sources. 
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